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We use a momentum-space hole-burning technique implemented via stimulated Raman transitions 
to measure the momentum relaxation time for a gas of bosonic atoms trapped in an optical lattice. 
By changing the lattice potential depth, we observe a smooth crossover between relaxation times 
larger and smaller than the bandwidth. The latter condition violates the Mott-Ioffe-Regal bound 
and indicates a breakdown of the quasi-particle picture. We produce a simple kinetic model that 
quantitatively predicts these relaxation times. Finally, we introduce a cooling technique based upon 
our hole-burning technique. 


Numerous quantum many-particle systems, such 
as electrons in metals, are surprisingly well modeled 
by gases of weakly interacting quasiparticles. This 
remarkable observation provides one of the most 
powerful paradigms in quantum many-body physics 
and underpins many calculational techniques. A 
quasiparticle is typically viewed as a composite 
formed from an elementary particle (e.g., an elec¬ 
tron) and a cloud of excitations (e.g., phonons or 
particle-hole pairs) [T]. The aggregate object be¬ 
haves like a particle: it has a well-defined quasimo¬ 
mentum k , energy e&, and mean-free path A&. The 
latter is the average distance the quasiparticle trav¬ 
els before being scattered into another direction. In 
1960, Ioffe and Regel pointed out that this picture is 
sensible only if the mean-free path is larger than the 
lattice spacing in the material [2], an idea which was 
further developed by Mott in 1972 [3]. In solids, esti¬ 
mates of the mean-free path from resistivity indicate 
that most metals satisfy this Mott-Ioffe-Regel (MIR) 
bound. Those materials that violate the MIR limit 
typically have other unusual properties [4 , 5 . For 
example, many have a resistivity that grows linearly 
with temperature T, in contrast to the T 2 prediction 
of Fermi-liquid theory. An equivalent expression of 
the MIR bound is that, within a quasiparticle pic¬ 
ture, the quasi-momentum relaxation time r asso¬ 
ciated with A*; cannot be shorter than h/W , where 
W is the bandwidth (i.e., the difference between the 
largest and smallest e^), and h = 2ivh is Planck’s 
constant. 

As we demonstrate, the disappearance of quasi¬ 
particles need not be a sudden event. For bosonic 
atoms trapped in an optical lattice, we find no sharp 
transition in the quasi-momentum relaxation time, 
which gradually decreases as the inter-particle inter¬ 
action strength is increased. For weak interactions 
(compared with tunneling) the relaxation is slow, in¬ 


dicating that the system is well described by quasi¬ 
particles. In contrast, for strong interactions, relax¬ 
ation is sufficiently fast to violate the MIR bound 
and suggest that the quasiparticle concept cannot be 
applied. These regimes are continuously connected, 
and we cannot demarcate an abrupt boundary be¬ 
tween them. This observation is consistent with the 
idea that the MIR bound is not a sharp inequality, 
but rather a scale describing a crossover between be¬ 
haviors. 

We study a gas of 87 Rb atoms trapped in a three- 
dimensional cubic optical lattice in the superfluid 
regime of the Bose-Hubbard model. The atoms are 
cooled to the lowest energy band of the lattice in the 
\l) = \F = 1 ,rriF = 0) hyperfine state. Significant 
occupation of high quasimomentum states is present, 
because the temperature of the gas is comparable to 
the bandwidth W = 12£, where t is the tunneling 
energy. As illustrated in Fig. 1(a), we use stimu¬ 
lated Raman transitions to drive the atoms to the 
first-excited band and |j") = \F = 2, mp = 0) state. 
Because the Raman wavevector difference A k = 
1.457r/d lies along the ^-direction of the lattice, the 
Raman transitions can only change the quasimomen¬ 
tum component q z by A q z = 2qs — hAk ~ 0.55 qs, 
where qs = hn/d and d = 406 nm is the lattice spac¬ 
ing. The strong dispersion of the first-excited band 
enables us to tune the frequency difference A uj be¬ 
tween the Raman beams so that atoms with q z « 0 
are resonant [6'. Fig. 1(b) shows images of quasimo¬ 
mentum distributions illustrating that these atoms 
are transferred to the |t) state. We subsequently 
eject these atoms from the trap using a 50 /is pulse 
of light resonant with the 5 Si/2, F = 2 —>> 5P 3 / 2 
transition. 

After removing the low q z atoms, we allow the 
remaining atoms to evolve in the lattice potential 
for a time thoid = (0 — 10) ms. We measure the 
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FIG. 1. Quasimomentum-selective stimulated Raman 
transition and procedure for measuring quasimomentum 
relaxation, (a) Energy e of 87 Rb atoms confined in a 
lattice potential. Only the ground band (n ass 0) in 
14-} and first-excited band (n = 1) in |j^) are shown. 
The Raman beams (red arrows), with Rabi rates 
and O 2 and detuned A from the transition to the 5 Pi/ 2 
excited electronic state, are used to selectively transfer 
atoms between bands, (b) Absorption images taken af¬ 
ter bandmapping and 20 ms of TOF. The images are 
obtained before (top) and immediately after (middle) 
the Raman pulse and subsequent to the resonant-light 
pulse (bottom). The reciprocal lattice vectors are indi¬ 
cated by k x ,k y , and k z . The projection of the FBZ onto 
the imaging plane is displayed with dashed lines. Atoms 
in the first-excited band appear outside the FBZ after 
bandmapping. 


quasi-momentum distribution after 0 id using band¬ 
mapping and time-of-flight (TOF) imaging |7J[8]. As 
shown in the inset of Fig. 2, the quasi-momentum 
profile rapidly relaxes after the q z ~ 0 atoms are 
removed. We use the mean squared residual of 
a fit to such images to determine the deviation 
from equilibrium at each thoid and measure the re¬ 
laxation timescale. The images are fit to a semi- 
classical model that describes the equilibrium quasi¬ 
momentum distribution of a non-interacting bosonic 
gas trapped in a combined lattice-parabolic poten¬ 
tial. The mean squared residual of the fit is r 2 = 
(OBij - rtijf / JAj 1, where OD ^ is the mea¬ 
sured optical depth, is the fit function [HUH, and 
the summations are over the indices i and j within 
a mask defined by the first Brillouin zone (FBZ). 
The residual r 2 is measured as thoid is varied at four 
different lattice depths: s = 4,5,6 and 8 Pr, where 
Er = h 2 n 2 /2md is the recoil energy, and m is the 
atomic mass. We observe that the quasimomentum 
relaxation is exponential in time, as shown in Fig. [2] 
for the particular case of s = 6Er. The measured 


relaxation time constant r determined from a sim¬ 
ple exponential decay for the four lattice depths is 
shown in Fig. it- 



FIG. 2. Sample relaxation data. The mean squared 
residual of the quasimomentum-distribution fit r 2 ver¬ 
sus hold time in the lattice Gold is shown for s = 6Er. 
Each data point corresponds to a single measurement, 
and the vertical axis of the plot has been re-scaled. The 
offset from r 2 — 0 is consistent with the residual from 
equilibrium images and is likely due to the failure of 
bandmapping at the edge of the FBZ (marked by the 
dashed white lines) and imaging noise. The red line is a 
fit to a single exponential decay used to determine the 
relaxation time constant r. 


The relaxation is comparable to or faster than the 
tunneling time h/t = (0.5 — 1.5) ms and similar to 
the Hubbard interaction time h/U = (0.4 — 0.2) ms 
for s = 6-8 Er , where U is the interaction energy 
for two atoms to occupy a site. Furthermore, r is ap¬ 
proximately two orders of magnitude faster than the 
time constant observed in the absence of the lattice 
(9]. The relaxation process speeds up for a stronger 
lattice potential, ultimately reaching just 200 /i s at 
5 = 8 Er. 

We compare the measured time constant to a 
short-range, two-body scattering calculation based 
on Fermi’s golden rule (FGR), treating the Hub¬ 
bard interaction term as a perturbation to the single¬ 
particle tight-binding Hamiltonian. We consider a 
thermal gas at equilibrium in a cubic lattice poten¬ 
tial and calculate the rate of scattering into the n^ =0 
state. We identify the relaxation time r with the in¬ 
verse of this rate, 
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where (n) = N(muj 2 d 2 / ^nksT) 3 / 2 is the filling- 
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FIG. 3. Quasimomentum relaxation measurements at different lattice potential depths, (a) The measured relaxation 
time r is shown for varied lattice potential depth s. Each point is determined using a fit to data such as those shown in 
Fig. 2, and the error bar displays the fit uncertainty, (b) Normalized relaxation rate h/tr versus U 2 /t 2 . The dashed 
line represents the scaling law predicted by Eq. 0 - The error bars include the fit uncertainty used to determine r 
and the standard deviation in measurements of T and N. The MIR bound is violated in the region marked in gray, 
which corresponds to h/rW > 1. 


weighted average lattice filling, ks is Boltzmann’s 
constant, N is the atom number, T is the tempera¬ 
ture at equilibrium, and Fit/ksT) ~ 0.3 is a numer¬ 
ical factor that does not depend strongly on t/ksT 
in the regime we probe 0- Fig HP compares the 
prediction from Eq. 0 with the measured relax¬ 
ation times. We show the normalized relaxation rate 
h/tr versus U 2 /t 2 with (n)=0.13, which is the aver¬ 
age value across all lattice depths determined by the 
lattice filling before relaxation. No free parameters 
were used for the theoretically predicted r, which is 
constrained by the known experimental values. 

Aside from an apparent saturation at small U/t, 
theory and experiment quantitatively agree. The 
source of the disagreement at low s is unknown and 
may arise from other relaxation processes, such as 
heating from the lattice light. Surprisingly, Eq 0 
appears to be valid even when U/t is as large as 7. 
Despite the fact that the theory is predicated on a 
quasiparticle picture, we find h/Wr ~ 4 at the high¬ 
est U/t, which is a violation of the MIR bound. We 
emphasize that there is no law preventing a “viola¬ 
tion” of the MIR limit—rather, the violation simply 
indicates that the elementary excitations are over¬ 
damped, and it is not sensible to describe them as 
quasiparticles. It is remarkable that the perturba¬ 
tion theory result can be so accurately extrapolated 
into this regime. 


The extraordinarily fast quasimomentum relax¬ 
ation we observe is key to enabling in-lattice cooling, 
which is an outstanding challenge for experiments 
focused on simulating models of strongly corre¬ 
lated electronic solids m • Several cooling schemes 
in optical lattices have been proposed 6, fTTHT?] . 
but experimental demonstrations have been limited 
[11 T9] . Notably, the quasimomentum degree of 
freedom has not been cooled directly and remains 
hotter than the Neel temperature in experiments 
with fermionic atoms m • To demonstrate the possi¬ 
bility of in-lattice quasimomentum cooling, we per¬ 
form a proof-of-principle experiment using Raman 
transitions to remove the most energetic atoms from 
the ground band, allowing the remaining atoms to 
equilibrate to a lower temperature. This technique 
is a momentum-space analog of evaporative cooling 
[21]. We repeat the process, thereby showing that 
this approach can be used as part of an iterative 
cooling method. 

The gas is initially prepared at the same temper¬ 
ature and atom number as for the relaxation experi¬ 
ment in a s = 4 Er lattice. For the cooling sequence 
(depicted in Fig. [4|, the Raman beams and reso¬ 
nant light are subsequently pulsed on for 400 /is and 
50 fi s, respectively. Each Raman pulse is designed to 
excite only atoms with high quasimomentum, while 
the condensate is left largely unaffected. Given the 
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spatial configuration of the Raman beams in our ex¬ 
periment, atoms are removed from one side of the 
FBZ in the k z direction, as shown in Fig. [4j At the 
end of a single cycle, a 1 ms delay is included to al¬ 
low the quasimomentum distribution to equilibrate. 
After all the cycles have completed, the gas is held 
in the lattice potential for 4 ms before bandmapping 
and imaging. 


Lattice 


Raman | 
Resonant light | 



lattice turn on 1 st cycle 2 nd cycle equilibration 




FIG. 4. Cooling of the quasimomentum distribution. 
The experimental sequence for cooling is shown schemat¬ 
ically at the top; the timeline is not to scale. The TOF 
images show the quasimomentum distribution immedi¬ 
ately before (i) and after (ii) atoms have been removed 
from the thermal component, (a) Quasimomentum pro¬ 
file along k z averaged over 3 — 4 images before (black) 
and after (gray-shaded) two cooling cycles, (b) Atom 
number in the condensate (a) and in the thermal com¬ 
ponent (#) as the cooling cycles are performed. (c),(d) 
The condensate fraction and the fit parameter (/3t) _1 
(which is monotonically related to the temperature) as a 
function of the cooling cycles. Plots (b)-(d) are obtained 
from fitting the TOF images to a semi-classical model. 
The error bars represent the standard error of the mean 
of the measurements averaged for each point. 

Figure^ shows the quasimomentum profile of the 
gas along k z before cooling and after performing 
two cooling cycles. We observe that the width of 
the thermal component shrinks and the condensate 
number grows as the cooling sequence is performed. 
For a quantitative analysis, we fit the TOF images to 
the non-interacting, semi-classical model used in the 
relaxation analysis, with an additional independent 


Thomas-Fermi profile for the condensate. We use 
the fit to determine the number of thermal and con¬ 
densate atoms (and the condensate fraction Nq/N ) 
and pt = t/ksT. 

The results of these fits are shown in Figs [4]3,c, 
and d. Figure [4 ]d shows that the thermal num¬ 
ber decreases as atoms are expelled from the trap, 
while the condensate number increases. The redis¬ 
tribution of atoms from the thermal component to 
the condensate is a sign of rethermalization during 
the cooling cycle. Moreover, the steady increase in 
the condensate fraction evident in Fig. 4c demon¬ 
strates that the entropy per particle is reduced dur¬ 
ing cooling, and this technique may therefore be used 
to reach new quantum phases that exist at lower 
entropy. This conclusion is further reinforced by 
Fig. [4]i , which shows a monotonic decrease in the 
temperature inferred from the fits to TOF images. 
Because bandmapping fails at the edge of the FBZ, 
a systematic error is made in the fitted f3t. Never¬ 
theless, the fitted (/?£)“ 1 is monotonically related to 
T [8], and the decrease in (/3£) _1 signals a reduction 
in temperature. 

A measure of the efficiency for any evaporative 
cooling scheme is a = dlogN/dlogT. A smaller 
value of a indicates more efficient cooling—fewer 
atoms are removed for the same change in T. For 
our method, a = 1.75 ± 0.04 based on a fit to the 
data shown in Fig. [4j This performance compares 
favorably with recent results for non-lattice gases, 
including a « 1.5 and 1.9 for “tilt” evaporation in 
a hybrid magnetic-optical trap [22j 23 and a ~ 2.7 
for dipole-trap evaporation of 87 Rb [24] . 

The ultimate limit to the lowest temperature 
achievable by any cooling method is determined by 
competition between cooling and heating rates— 
cooling ceases when the two are equal [10]. The 
heating rate in optical lattices is primarily deter¬ 
mined by momentum diffusion resulting from the 
interaction between the light and atoms [25]. This 
effect can be minimized by detuning the light far 
from any electronic transition. In our experiment, 
the heating rate induced by the lattice light is 
0.15 E r /s=25 pK k R /ms at s = 4 E R , while the 
cooling power based on the data shown in Fig. §i 
is approximately 9 nK k R /ms (corresponding to 0.6 
t/ ms). This extraordinary cooling power is possi¬ 
ble because of the high thermalization rate. In the 
regime we explore, heating from the lattice is there¬ 
fore not a limitation to the cooling method. 

In principle, more cooling cycles could have been 
carried out in our experiment. However, drift in the 
initial temperature of the gas limited our ability to 
optimize the pulse sequence. The cooling efficiency 
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can be improved by addressing both sides of the FBZ 
simultaneously using a small Raman wavevector A k 
compared with ir/d. Moreover, further enhancement 
in efficiency can be achieved by adding more pairs 
of Raman beams to target atoms in the transverse 
directions to k z j9]. The ultimate cooling limit for 
this technique, as applied to bosonic atoms, is set by 
the selectivity in quasimomentum together with the 
finite momentum spread of the condensate. 

In conclusion, we have developed a novel tech¬ 
nique for measuring quasimomentum relaxation 
times in an optical lattice. We discovered a viola¬ 
tion of the MIR bound and showed how our method 
can be adapted to cool any atomic species. Further 
studies of relaxation in this system may contribute 
to our understanding of transport in materials that 
also display a violation of the MIR bound [4] and 
thermalization in closed quantum systems (26[. In 
the future, the cooling technique could be applied 
to fermionic atoms trapped in optical lattices in or¬ 
der to reach exotic quantum states that may exist 
at low entropy per particle. 
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